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Abstract 

We study CP-violating effects in — > ^l^Kg decay witliin minimal 
supersymmetric models with spontaneous CP-violation. We find that 
the CP-asymmetry predicted by the Standard Model in this decay, 
sin 2/? > 0.4, cannot be accommodated in these models without vio- 
lating the bound on the neutron electric dipole moment. This result 
holds for NMSSM-like models with an arbitrary number of sterile su- 
perfields. Further implications of the scenario are discussed. 



1 Introduction 



The origin of CP-violation is one of the most profound problems in parti- 
cle physics. In the Standard Model, all observable CP-violating effects in 
the kaon system can be successfully explained via the Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism [1]. However, the physical principles lying be- 
hind CP-violation are still not understood. 

One of the more elegant approaches to the problem of CP-violation is 
based on the possibility of spontaneous T-breaking in multi Higgs doublet 
systems [2]. Supersymmetric models can provide such systems and thus are 
a natural setting to implement this idea. Spontaneous CP- violation (SCPV) 
in susy models has drawn considerable attention [5-10] due to the following 
attractive features: 

• CP-phases become dynamical variables, 

• CP-symmetry is restored at high energies, 

• it allows to avoid excessive CP-violation inherent in susy models. 

In this letter, we consider the implications of SCPV in minimal susy models 
for CP-asymmetries in 5 — > ipKg decays. It is well known that the Standard 
Model predicts large CP- violation in these decays, namely sin 2/5 > 0.4 where 
f3 is the one of the angles of the unitarity triangle [21]. Currently this CP- 
asymmetry is being studied experimentally by the CDF collaboration. Even 
though the statistics does not allow to make definite statements about the 
validity of the SM predictions at the moment, large CP- violation in this 
decay has been hinted. The purpose of this paper is to determine whether 
a large CP-asymmetry, sin 2/3 > 0.4, can be explained in susy models with 
spontaneously broken CP. 

2 CP- Asymmetry in B ^ ipKs Decay and 
Spontaneous CP- Violation 

In this section we will consider minimal susy models with spontaneously 
broken CP and obtain the lower bound on the CP- violating phase as dictated 
by sin 2f3 > 0.4 . 

It has been shown that the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) is the simplest susy model which allows spontaneous CP- 
violation while being consistent with the experimental bound on the lightest 
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Higgs mass [8]. In the most general version of the NMSSM with the super- 
potential 

W = XNH,H2 - ^iV^ -rN + fiH.H^ + W^/™„ , (1) 

SCPV can occur already at the tree level thereby avoiding the Georgi-Pais 
theorem [6]. Note that even though SCPV in the MSSM is allowed theoret- 
ically [5], such a scenario is ruled out by the LEP constraints on the axion 
mass [7]. 

We will assume that all CP-violating effects result from the complex Higgs 
VEV's 

{H',)=vu {H',)=v2e'',{N)=ne'^ . (2) 

The relevant interactions^ for one generation of fermions (after spontaneous 
SU{2) X f/(l) symmetry breaking) can be written as follows [9,12]: 

£ = -KH^URUL-hdH^dRdL-gW^PLdul + hddPLH^UL 

+ huH^'PLd u*p. + h.c. , (3) 

Crm. = -g{viHPLW + V2e-'^WPLH) + e-'''-Km^;:l\*j,UL 

+ e-^^'^h^mfj^drjL + h.c. , (4) 

where hu^ denote the Yukawa couplings and Ku,d are certain functions of 
the Higgs VEV phases p and ^. In what follows, we assume, for the sake 
of definiteness, that = \kii\ = \k\, m^^j^ = m^'^^ = m'j^j^, and that 
K and mj^pi^ are generation independent. Equation (4) represents the wino- 
higgsino and left-right squark mixings, with the former being responsible for 
the formation of the mass eigenstates - charginos. However, following Ref. 
[9], we will prefer the "weak" (wino-higgsino) basis over the mass (chargino) 
one. 

Information about the angle (3 of the unitarity triangle was extracted from 
the decay rate evolution 

r(S°[5°](t) ^ iIjKs) oc e-^'(l - [+] sin 2/? sin Amt^ , 
^The complete list of interactions can be found in [12]. 
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where Am is the Bl — Bh mass difference. On the other hand, the angles of 
the unitarity triangle can be expressed in terms of the mixing (0m) and decay 
(0£)) phases which enter the B — B mixing and h — > qqQ decay diagrams. 
For the process under consideration, 

sin 2(3 = sin(20D + 0m) • (5) 

This relation is theoretically clean since it does not involve hadronic uncer- 
tainties and can serve as a sensitive probe for physics beyond the Standard 
Model. At the present time, the CKM entries are not known precisely enough 
to make a definite prediction for (3. However, it is known that sin 2(3 must 
fall between 0.4 and 0.9 in order for the CKM model to be consistent [21]. In 
our model, we will impose this condition together with Eq. (5) to obtain the 
lower bound on the CP-violating phases appearing in the decay and mixing 
diagrams. 

Let us now proceed to calculating the CP- violating effects in B — B mixing. 
Figure 1 displays what we believe to be the most important contributions 
to the real part oi B — B mixing. These include the Standard Model box 
and wino superbox diagrams. It has been argued [9] that all significant 
CP-violating effects result from complex phases in the propagators of the 
superparticles. For the K — K system, one loop diagrams involving complex 
phases in the wino-higgsino and left-right squark mixings can lead to the 
correct values of e and e' [9]. The analogous B — B diagrams are shown in 
Figure 2 (in the case oi K — K mixing, the diagram in Fig. 2b is super-CKM 
suppressed and can be neglected). Besides the above mentioned contribu- 
tions, there is a number of other contributions to B — B mixing which can 
be classified as follows: 

1. Higgs boxes, 

2. gluino superboxes, 

3. neutralino superboxes. 

The gluino contribution can be neglected since the gluino is likely to be very 
heavy: rrig > 310 GeV [13]. Further, the neutrahno analogs of Figs.2a,b have 
to involve at least two powers of hb and, thus, are suppressed by {rrih/mw)'^- 
The same argument applies to the charged Higgs contribution [10]. On the 
other hand, the neutralino and Higgs contributions to Re {B — B) can be 
significant. However, they interfere with the SM contribution constructively 
[14] and can only reduce the B — B mixing weak phase. Since our purpose 
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is to determine the lower bound on this phase as dictated by sin 2/3 > 0.4, 
we can safely omit these corrections. 

Therefore, for our purposes it is sufficient to retain the SM and chargino 
contributions only. Furthermore, note that we may concentrate on the {V — 
A) X (V—A) interaction solely since 4-fermion chargino-generated interactions 
involving the right handed quarks arc suppressed by powers of miy/niw. As a 
consequence, hadronic uncertainties and QCD corrections will not affect our 
results since those will cancel in the expression for the phase. The resulting 
4-fermion interaction can be expressed as 

Oab=2 = {ksM + ksusY + e''l2 + ze-''l2 + e^''h + z''e-^''U) 

X dYPLbd-f^Pib , (6) 

where kgM and ksusY arc real couplings induced by the diagrams in Fig. la 
and Fig. lb, respectively. The CP- violating couplings e*'^/2 and 6^**^/4 result 
from the diagrams with two and four complex mixings shown in Figs. 2a and 
2b, respectively (^2,4 are defined to be real). It is important to note that along 
with the diagrams explicitly shown in Fig. 2, there are also "cross" diagrams 
in which the positions ofii and are interchanged. Such graphs contribute 
with opposite phase and may seem to lead to a complete cancellation of the 
imaginary part of the coupling. However, this cancellation is only partial [9] 
owing to the fact that the higgsino vertex is, generally speaking, different 
from that of gaugino. Such partial cancellation is accounted for by a variable 
factor z {0< z<l). 

The Standard Model contribution is given by [15]^ 

ksM = ^mlH{xt) {VM' , (7) 

with H{x) being the Inami-Lim function. To estimate the superbox con- 
tributions, we may treat the gaugino and higgsino as particles of mass 
with (perturbative) mixing given by Eq.(4). In this approximation, 
the chargino propagator with a mixing insertion (Fig. 2) is proportional to 
g(vi + f 26"^^) (1,21^^ )2 ■ The resulting 4-fermion couphngs are given by 



ksuSY = 9 ~^ (KfcKd)^ 



1287r2 m? ^ {y - 1) 



y + 1 -Iny 

1/ - 1 



(8) 



We do not show the QCD correction factor exphcitly. 
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,(10) 



where y = rn^/m'^ ; and m^y denote the top squark and the chargino 
masses, respectively, and V is the squark analog of the CKM matrix. In 
these considerations, the top squark contribution is believed to play the most 
important role. The influence of the c- and w-squarks is taken into account 
through a variable super-GIM cancellation factor ( {0 < ( < 1). Such a 
factor is associated with every squark line on which summation over all the 
up-squarks takes place. Since masses of the top and c, w-squarks are expected 
to be very different due to the large top Yukawa coupling (as motivated by 
SUGRA), the natural value of ( would be of order unity {( is deflned to 
vanish in the limit of degenerate squarks). 

To derive the lower bound on the phase 6, we may replace Vi + V2e~^P in 
Eqs.(9) and (10) by its "maximal" value \/2ve~'^'^ with v deflned in the usual 



way: f = yuf + ~ lYAGeV. Apparently, 5 < |p| + |k| . 

Let us now determine the weak phases (pM and (po. It follows from Eq. (6) 
that 

^ ^ l2{l — z) smS + U{l — z'^) sin2(5 

ksM + ksusY + ^2(1 + z) cos5 + ^4(1 + z"^) cos 25 ' 

On the other hand, the decay phase (pD is negligibly small [10]. Indeed, in 
our model, the only source of CP-violation in the process b — > ccs is the 
superpenguin diagram with the charginos and squarks in the loop|^. How- 
ever, this diagram is greatly suppressed as compared to its CP-conserving 
counterpart, W mediated tree level decay, due to the loop factors and heavy 
squark propagators. Also, unlike for the kaon decays, there is no A/ = 1/2 

■^Another possible contributor, Higgs- mediated tree level decay, is suppressed by the 
quark masses. 
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enhancement of the {V — A) x {V + A) interactions. As a result, the weak 
decay phases can be neglected. This also means that direct CP-violation in 
our model is neghgibly small as compared to that in the Standard Model 
(unless there is no tree level decay mode). 
According to Eq.(5), the angle /? is determined by 

sin 2/3 = sin 0m ■ (12) 

Then the experimental bound sin 2/3 > 0.4 can be translated into 

tan (/)M > 0.44 . (13) 

This, in turn, leads to a lower bound on the phase S which can be obtained 
numerically from Eq.(ll). Note that Eq.(ll) is free of hadronic uncertanties 
and QCD radiative corrections. 

3 Numerical Results 

In this section we will discuss implications of Eq.(13) and its compatibility 
with the upper bound on the NEDM. 

It is well known that the tight experimental bound on the NEDM imposes 
stringent constraints on the CP-violating phases which appear in extensions 
of the Standard Model. In our model, the largest contributions to the NEDM 
are shown in Fig. 3. Barring accidental cancellations, one can constrain the 
CP-phases entering the higgsino-gaugino and squark left-right mixings via the 
chargino (Fig. 3b) and gluino (Fig. 3a) contributions to the NEDM, respec- 
tively. Consequently, the phase S appearing in the B — B mixing becomes 
bounded due to 

sin 5 < I sin k,\ + \ sin p\ . (14) 

This requires d to be of order 10~^-10~^ [9,16]. 

On the other hand, Eq.(13) imposes a lower bound on 6. As seen from 
Eqs.(8)-(10), this bound depends strongly on the super-CKM matrix. We 
consider the following possibihties: 

1. the super-CKM matrix duphcates the CKM one, Vtd ~ Vtd; 

2. the super-CKM mixing is enhanced, Vtd ~ Vtd/ sin^^c*; 

3. the super-CKM mixing is doubly enhanced, Vtd ~ Vj^/ (sin 6'c)^- 
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In all of these cases we assume Vf?, ~ 1. The first possibility implies that the 
super symmetric contribution to B — B mixing is suppressed as compared to 
the CP-conserving Standard Model box diagram. As a result, we find that 
the constraint tan0jv/ > 0.44 cannot be satisfied for any 6 even assuming 
light (100 GeV) squarks. For the same reason, we are bound to consider 
only light (100 GeV) chargino and maximal left-right squark mixing, ttllji = 
rriq. On the other hand, the third option is unrealistic since it leads to an 
unacceptably large stop contribution to the Ks — Kl mass difference unless 
the top squark mass is around 1 TeV. Therefore, we are left with the second 
possibility which we will examine in detail. From now on we assume that 
Vtd — Vtd/ sin ^(7, rriy^r — lOOGeV, miR — rriq, and will study the behavior 
of the lower bound on 5 as a function of the remaining free parameters - 
z, (, tan P, and nig. 

Fig. 4 displays a typical picture showing inconsistency of the modelQ The 
lower bound exceeds the upper bound by one or two orders of magnitude. 
Moreover, the region allowed by the CP-asymmetry in B "ipKs is re- 
stricted to the left upper corner of the plot. The reason for that can be 
easily understood. Indeed, if the squarks are heavy, the magnitude of the 
susy contribution is negligible as compared to that of the CP-conserving SM 
box and sufficient CP- violation cannot be produced regardless of how large 
the phases are. 

Let us now consider the effect of each of the variable parameters.^ 

1. — dependence {Fig. 5). 

The lower bound on 6 relaxes as we introduce the super-GIM cancellation. 
This occurs due to the increasing share of the CP-violating diagram in Fig.2b. 
However, theoretically one expects ( to be of order unity due a large differ- 
ence between the stop and other squarks masses. 

2. tan/? — dependence {Fig.6). 

The gap between the lower and upper bounds widens drastically as tan/? 
increases. Recalling that hu = . ^ and hd = ^ ^^'^ — it is easy to see 

^The displayed bound on the NEDM was calculated for the gluino mass in the range 
300 — 500 GeV using the standard formulas [16]. 

^ We are assuming that 5 belongs to the second quarter. Even stricter lower bounds 
on sin (5 can be obtained for 5 in the second quarter due to a partial cancellation in the 
numerator of Eq.(ll). 
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that for large tan/? the NEDM constraint becomes stricter due to the large 
hd whereas the CP- violating contributions to B — B mixing, proportional to 
powers of /it, decrease. We do not consider the case tan/3 < 1 because of 
the SUGRA constraints and the breakdown of perturbation theory in this 
region. 

3. z — dependence {Fig. 7). 

Apparently, the incorporation of a partial cancellation {z > 0) among the 
CP-violationg contributions makes the lower bound rise. One expects the 
natural value for z to be around 1/2. 

In all of these cases the regions allowed by the NEDM and B — > ipKs 
are at least an order of magnitude apart.Q Furthermore, heavy squarks (> 
400 GeV) are prohibited by large CP-asymmetries observed in i? ^ 
This condition is quite restrictive and may alone be sufficient to rule out the 
model in the near future (even if large CP-phases were allowed). 

It should be mentioned that, in the limit of large tan/?, the CP-violating 
neutralino and Higgs contributions to B — B mixing become more impor- 
tant. However, the CP-phases are severely constrained in this case (~ 10~^). 
Therefore, these contributions do not lead to any considerable modifications 
of our analysis. 

4 Further Discussion 

The model under consideration has further implications for B-physics. For 
instance, in this model, direct CP-violation is greatly suppressed in all tree- 
level allowed processes as compared to what one expects in the Standard 
Model. This provides another signature testable in the near future. 

The other angles of the unitarity triangle, a and 7, can be determined in 
a similar manner from, for example, B^ — ^ tt+tt" and Bg — >■ pKg decays [11]. 
However, in the Standard Model, one cannot determine these angles from 
the decay rate evolution and relations analogous to (5) precisely enough due 
to considerable penguin contributions. To eliminate their infiuence, one can 
use isospin and SU{3) relations [18]. For instance, in order to determine a, 
one needs to know the rates of B^ —>■ tt+tt", B^ — > vr^vr^, and B^ —>■ 7r+7r° 

^Note also that small CP-phases are disfavored by the bound on the lightest Higgs 
mass [8]. 
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along with the rates of their charge conjugated processes. Then a can be 
found from certain triangle relations among the corresponding amplitudes. 
In our case, however, this analysis becomes trivial due to a vanishingly small 
interference between the tree and superpenguin contributions. As a result, 
a can be read off directly from the analog of Eq.(12). The angles extracted 
in such a way normally do not exceed a few degrees (modulo 180") and do 
not form a triangle [10]. 

In our model, the angles of the unitarity triangle do not have a process- 
independent meaning. Indeed, contrary to the Standard Model, they cannot 
be extracted from direct CP-violating processes simply because such pro- 
cesses are prohibited. Moreover, these angles are not related to the sides of 
the unitarity triangle. 

Another important consequence of the model is that the CKM matrix is 
real and orthogonal. This, of course, is also true for general (nonsupersym- 
metric) two Higgs doublet models with FCNC constraints, in which CP is 
broken spontaneously [19]. As a result, and | sin6'cVcfe| must form 

a flat triangle. Presently, such a triangle is experimentally allowed provided 
new physics contributes significantly to Ams^ [20]. To determine the status 
of these models, a more precise determination of |V^b| and \ytd\ is neces- 
sary. Orthogonality relations among other CKM entries are less suitable for 
probing this class of models due to the small CKM-phases involved. 

We have considered in detail the case of the NMSSM. The results, however, 
remain valid for NMSSM-like models with an arbitrary number of sterile su- 
perfields iV. Indeed, since the iV's do not interact with matter fields directly, 
an introduction of a sterile superfield does not affect the way CP-phases en- 
ter the observables. The CP-violating effects will still be described by the 
diagrams in Figs. 2 and 3. Therefore, the argument we used also applies to 
this more general situation: the CP-phase in left-right squark mixing can 
be constrained via the gluino contribution to the NEDM, whereas the phase 
in the gaugino-higgsino mixing can be constrained via that of the chargino 
(assuming no accidental cancellation). In the same way, we find that the 
lower and upper bounds on the CP-phases are incompatible. A nontrivial 
extension of the model, in which a cancellation among various contributions 
to the NEDM can be well motivated, is necessary to rectify this problem. [| 

^ The possibility of such a cancellation in certain susy models has recently been con- 
sidered by a few authors. It is, however, unclear whether this cancellation can be made 
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5 Conclusions 



We have analyzed the CP-asymmetry in S — > ipKg decay within minimal 
supersymmetric models with spontaneous CP-violation. We have found that 
the CP-asymmetry required by sin 2/9 > 0.4 can be accommodated in these 
models only if the following conditions are met: 

1. left-right squark mixing is maximal, m^i? ~ rrig, 

2. super-CKM mixing is enhanced, Vtd — Vtd/ sinOc, 

3. the chargino is relatively light, ~ 100 GeV, 

4. the t-squark is lighter than 350 — 400 GeV. 

Even if this is the case, the required CP-violating phases are larger than 
those allowed by the bound on the NEDM by one or two orders of magnitude. 

We conclude that the model under consideration cannot accommodate a 
large CP-asymmetry in i? — > i/jKg while complying with the bound on the 
NEDM. Thus, if the Standard Model prediction gets confirmed, the model 
will be recognized unrealistic. The result holds true for models with an 
arbitrary number of sterile superfields. To reconcile theory with experiment 
one would have to resort to essentially nonminimal scenarios in which large 
CP-violating phases are naturally allowed. 

Note also that since, in this approach, CP-violation is a purely supersym- 
metric effect, the model requires the existence of a relatively light t-squark. 
This may conflict with the Tevatron constraints on supersymmetry in the 
near future (see, for example, [13]). 

We have also discussed other testable predictions of the model. Among 
them are the suppression of direct CP-violation and orthogonality of the 
CKM matrix. 

The author is grateful to L. N. Chang and T. Takeuchi for discussions and 
reading of the manuscript, and to T. Falk for useful comments. 



natural (see [17] and references therein). 
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Figure Captions 

Fig. la,b Most important contributions to B — B mixing. 

Fig. 2a,b Dominant CP-violating contributions to B — B mixing (all pos- 
sible permutations are implied). 

Fig. 3a,b Contributions to the NEDM allowing to constrain the wino- 
higgsino and squark left-right mixing phases. 

Fig. 4 Regions allowed by the CP-asymmetry in i? — > iJjKs (upper) and 
the bound on the NEDM (lower). Typically, the region allowed hj B ^ 'ipKg 
is much smaller than shown here due to the partial cancellation (see Fig. 7, 
z^l/2). 

Fig. 5 Lower bound on sin 5 as a function of the super-GIM cancella- 
tion parameter 1 - ( = 1, 2 - C = 1/2, 3 - C = 1/4- 

Fig. 6 Lower bound on sin 5 as a function of tan/3: 1 - tan/3 = 1, 
2 - tan /3 = 25. 

Fig. 7 Lower bound on sin 5 as a function of the partial cancellation 
parameter z: 1 - z — 0, 2 - z — 1/4, 3 - z — 1/2. 
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